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We study the effect of gapless quasiparticles in a d-wave superconductor on the T — end 
point of the Kosterlitz-Thouless transition line in underdoped high-temperature superconductors. 
Starting from a lattice model that has gapless fermions coupled to 3D XY phase fluctuations of the 
superconducting order parameter, we propose a continuum field theory to describe the quantum 
phase transition between the d-wave superconductor and the spin-density-wave insulator. Without 
fermions the theory reduces to the standard Higgs scalar electrodynamics (HSE), which is known 
to have the critical point in the inverted XY universality class. Extending the renormalization 
group calculation for the HSE to include the coupling to fermions, we find that the qualitative effect 
of fermions is to increase the portion of the space of coupling constants where the transition is 
discontinuous. The critical exponents at the stable fixed point vary continuously with the number 
of fermion fields N, and we estimate the correlation length exponent (v = 0.65) and the vortex 
field anomalous dimension (77$ = —0.48) at the quantum critical point for the physical case N = 2. 
The stable critical point in the theory disappears for the number of Dirac fermions N > N c , with 
N c ~ 3.4 in our approximation. We discuss the relationship between the superconducting and 
the chiral (SDW) transitions, and point to some interesting parallels between our theory and the 
Thirring model. 



I. INTRODUCTION 

All high-temperature superconductors are well estab- 
lished to have d-wave symmetry of their order parameter 
0,0, and therefore to posses gapless fermionic excita- 
tions near the four points in the momentum space where 
the order parameter vanishes. Another distinguishing 
feature of high-temperature superconductors is that, to a 
very good approximation, they may be considered to be 
quasi-two-dimensional when underdoped || . This leads 
one to expect the superconducting transition at finite 
temperatures (T ^ 0) in the underdoped regime to be of 
the Kosterlitz-Thouless type j|. Indeed, there is mount- 
ing experimental support in favour of free vortices near 
and above the superconducting critical temperature ||, 
||. Building on these two observations, one of us has 
recently shown that in the minimal theory that couples 
vortex and quasiparticle degrees of freedom || the result 
of phase disordering in a d-wave superconductor (dSC) 
at T = is nothing but the antiferromagnetic (more pre- 
cisely, the incommensurate spin-density wave (SDW)) in- 
sulator, the phase that is believed to be neighbouring the 
dSC in the phase diagram of all high-temperature super- 
conductors. This result follows from realizing that inte- 
gration of the unbound vortex degrees of freedom yields 
the theory for the low-energy fermions that is closely re- 
lated to the two-component three dimensional quantum 
electrodynamics (QED 3 ) ||, @ in which the dynamical 
breaking || of an approximate " chiral" symmetry trans- 
lates into SDW ordering ]lC| j. 

One can show |TlJ that gapless quasiparticles do not 
change the Kosterlitz-Thouless universality class of the 
T =/= superconducting transition, as one could also ex- 
pect on quite general grounds. At T = 0, on the other 
hand, the effect of gapless fermionic excitations on the 
order parameter is often less benign, and may change 



the universality class, or even the order, of the phase 
transition In this paper we study the quantum 

(T = 0) phase transition in the theory of coupled gapless 
quasiparticles and vortex-loop excitations in the phase of 
the superconducting order parameter. On the insulating 
side, i. e. when the vortex-loops are unbound, our theory 
reduces to the effective QED 3 description of underdoped 
cuprates , [|| ■ Assuming that the dynamics of the su- 
perconducting phase at T = is governed by the three 
dimensional (3D) XY model, we show that the effect of 
gapless quasiparticles is then to introduce an additional 
massless gauge field into the Higgs scalar electrodynam- 
ics (HSE), which represents the dual theory for vortex- 
loop excitations. The effect of this extra gauge- field on 
the quantum critical point is controlled by the number 
of fermionic fields N, and for N = we recover the re- 
sults for standard HSE jl3| . Using the approximate fixed 
dimension renormalization group (RG) equations for the 
HSE, known to be qualitatively and even quantitatively 
reliable when N = we show that there exists a 
critical number of fermionic fields N c « 3.4, above which 
there are no longer stable fixed points in the theory, and 
the flow in the space of coupling constants is always to- 
wards the region with a negative quartic term coupling. 
We interpret these runaway flows as the fluctuation in- 
duced first-order transition. Similar conclusion is reached 
in a more controlled calculation in 4— e dimensions, where 
we find that the critical number of complex components 
of the scalar field above which a second order phase tran- 
sition can occur jf4| increases when fermions are present. 
This suggests that gapless quasiparticles tend to desta- 
bilize the inverted XY phase transition in favour of a 
first order one. For the physical case N = 2, however, 
we still find a critical point within our calculation, and 
we assume it controls the quantum phase transition for a 
strongly type-II HSE. Its critical exponents are estimated 
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to be v = 0.65 and 77$ = —0.48. 

The paper is organised as follows. In the next section 
we define the lattice theory one can use to model the 
system of gapless fermionic excitations coupled to vor- 
tex loops. In the sec. Ill we propose the continuum 
description of the lattice theory and discuss some of its 
features and possible problems. We then proceed to ex- 
tend the fixed dimension RG calculation of ref. onto 
the present problem in sec. IV. A supporting computa- 
tion in 4 — e dimensions is provided in sec. V. In sec. VI 
we discuss a simpler version of our theory. Conclusion 
and discussion of our results are presented in sec. VII. 
Technical details are given in Appendices. 

II. THE LATTICE MODEL 

We begin with the simplest lattice model that has gap- 
less fermions coupled to the XY variables: 

Z= f Y[(dUx)d Xl (x)d a (x)) (1) 

*' X 

I 

E exp-(5i[x,a] + S 2 [a,mi,m B ]), 

m A (x),m B (x) 

where 

N/2 

SiLval^E E ( 2 ) 

and 

S 2 [a,m A ,m s ] = E^[ Vx KW+^Wlf ( 3 ) 

X 

+ia(x) ■ [V x (m A (x) - m B (x))}}. 

Here x labels the sites of a three dimensional quadratic 
lattice, /i = 0,1,2, and Vx denotes the lattice curl. 
Prime on the sum in (1) means that V • m AtB = 0, 
where V is the lattice gradient, and m AB are integer 
vector variables. rj^{x) are the standard Kawamoto-Smit 
phases {m{x) = 1, m(x) = (-1) X1 , Vo(x) = (-l) ai+9i ») 
introduced to ensure relativistic covariance in the contin- 
uum limit Jl5| , [fl6f . Xi ( x ) are one component Grasmann 
variables. To be specific, we are using staggered fermions, 
and will, therefore, consider N to be even. 

Detailed derivation of the above partition function may 
be found in fjj. In the continuum limit it describes N 
four-component Dirac fermions which represent the 
neutral gapless spin-1/2 excitations (spinons) that can be 
identified in the d-wave superconducting state, coupled 
to XY degrees of freedom. Two species of integers (A and 
B) are the result of making the Franz- Tesanovic transfor- 
mation on the original electrons || , in which spin up and 
spin down electrons need to be transformed differently. 



Notice that in the absence of fermions (when N = 0) the 
integration over the gauge-field a enforces the constraint 
1x1,4(3:) = m B (x), so that the partition function reduces 
to the well-known current-loop representation of the XY 
model [^7). The parameter K is the XY phase stiffness, 
which may be understood as proportional to doping [fl8| , 
at least in the low doping regime of high-temperature 
superconductors. 

The problem we wish to address here is how the XY 
universality class is modified by the presence of gapless 
fermionic excitations in the d-wave superconductor. The 
case which is of particular interest is when N = 2, corre- 
sponding to four nodes of the d-wave order parameter in 
a single-layer superconductor. 

The lattice theory (1) can in principle be simulated on 
a computer, and the nature of the phase transition as the 
parameter K is varied studied that way. In this paper 
we will make some further transformations that will en- 
able us to suggest a continuum field theory which should 
have the phase transition in the same universality class 
as the lattice model. For that purpose we first rewrite 
the above partition function by using the standard dual 
(vortex loop) variables 1,4 and \ B as 

Z=lim / \\{d Xl {x)d Xl {x)da{x)dA A {x)dA B {x)) (4) 

y^Oj 11 
/ 

E exp-{S'i[x,a] + S 2 [a,A A ,A B ] + 

1a(x),\b(x) 

E i2ir(l A (x) ■ A A (x) + l B (x) ■ A B (x)) 

X 

+y(l\(x) + ll(x))}. 

When the "chemical potential" for the vortex loops y = 
the summation over integers l A (x) and l B (x) forces 
A A (x) and A B (x) to take strictly integer values, and 
one recovers the previous expression for the partition 
function. The approximation we will make is to take 
y small but finite. This approximation is often made 
in studies of the 2D |l9) and the 3D |2(J XY model, 
and is known not to change the universality class of the 
transition in that case |2lf . In present case the inno- 
cence of such an approximation is less clear. First, when 
y = the action is symmetric under the transformation 
a^(x) — ► a^(x) + 2irn fJ _(x) with n^(x) integer, and the 
summation over integers would in principle yield 

a compact gauge theory. This symmetry is broken when 
y ^ 0, and one therefore replaces in this way a compact 
gauge theory with a non-compact one. The difference be- 
tween the two, when the gauge-field is coupled to matter, 
is at present far from clear |22[] , and we will unfortunately 
have nothing further to add on this issue. Second, in the 
quenched limit N = the partition function for y ^= 
does not exactly reduce to the XY model any longer. 
Consequently, our continuum dual theory, which will as- 
sume y at N — 0, will not automatically be in the 
inverted XY universality class, as it should in the absence 
of fermions. We will, however, propose a way how to deal 
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with this problem in the next section. 

With these two caveats in mind we can proceed to- 
wards the continuum limit by observing that the last ex- 
pression is up to the Villain transformation |p3| equiva- 
lent to 

Z — lim / Y\(dxi(x)dxt{x)da.(x)dAA(x)dAB(x) (5) 

y^Oj 11 

d9 A {x)d9 B {x)) exp —{Si [Xt a ] + ^[a, A^, A.b] 
-— ^{cos(6»a(x) - ^(x + z>) - 2kA a , v {x)) 

+ cos{9 B {x) - 9 B {x + P)- 2nA B , v {x))}. 

Relaxing the y — > limit enables one to recognise the 
continuum limit of the lattice theory given in Eq. (5), 
which is what we do in the next section. 



III. DUAL CONTINUUM THEORY 

We write down a generalized continuum theory using 
n-component complex scalar fields. This model will be in 
the same universality class as the lattice {y ^ 0) model 
in Eq. (5) for n = 1, and bi — 62 ^ e 2 (strong type-II 
limit): 

C = Ch + Cf + C c , 

N 

C h = I (V - iA+ - zA_) *x| 2 + < |*i| 2 
+ |(V-iA + +zA_)* 2 | 2 +<ol*2| 2 
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c r = 



fl*ll 2 l*2l 2 

(V x A + ) 2 
ia • V x A_ 



(6) 



'J., are the four- component spinon fields |7| that represent 
the gapless quasiparticle excitations, and which are the 
continuum limit of the Xi lattice fermions. A and A + 
are the continuum versions of ^(A^ — As) and tt{Aa + 
Ab)- The two scalar fields 3>i and $ 2 are the "soft- 
spin" representation of the dual angles 8a and 9b in the 
continuum limit. There are two couplings that describe 
the interactions between 4>i and $ 2 fields, bi and bi. 
We have also introduced the parameter A to control the 
strength of the coupling to spinons. The physical case 
corresponds to A = 1. Finally, tuning the parameter K 
in the lattice theory should correspond to varying to in 
the Eq. (6). 

At the mean-field level the above continuum theory 
has a single phase transition from the insulating phase of 
unbound vortex loops (4>i) = (^2) = (3?) 7^ 0, when 
to < 0, to the superconducting phase (bound vortex 
loops) ($1) = (3>2) = 0, when to > 0. To understand 



the fate of fermionic degrees of freedom in both phases 
one needs to integrate out the gauge fields A_ and A + : 

1) In the superconducting phase A + is decoupled from 
a, and therefore decoupled from spinons. Integration 
over A_ gives a constraint V x a = 0, and therefore the 
spinons are free. This corresponds to sharp quasiparticle 
excitations in the dSC. 

2) In the insulating phase there is a Meissner effect for 
both A + and A , and the term ~ 2|(*)| 2 (A 2 h + A 2 _) 
appears in the action. A + is still decoupled from a, but 
the integration over A_ now yields 



N 



(V x a) 2 
32tt 2 |(*) 



(7) 



This is three dimensional quantum electrodynamics 
{QED 3 ) @, || . For N < N c ps 3 the "crural" symmetry 
|g| is spontaneously broken, and the "chiral condensate" 
(ipip) oc |(3>)| 2 appears. Going back to the electronic ori- 
gin of the QED3 one finds that the chiral condensate is 
precisely the SDW order parameter 0. 

The continuum model has three local U(l) gauge sym- 
metries, one for each of the three gauge fields. These 
symmetries should insure that the theory is renormahz- 
able, and that A_ and A+, for example, do not acquire 
a mass in the superconducting ((€>) = 0) phase. To pre- 
serve this gauge symmetry, however, one needs to keep 
two types of the fluctuating $ fields: if one would simply 
set $1 = €>2 hi Eq. (6) the model would lose the U{1) 
symmetry associated with the A_ field. This would lead 
to a theory that is not renormahzable; the A_ would 
becomes massive in the superconducting phase, as the 
Ward identities which usually prevent this from happen- 
ing would be no longer respected. 

In absence of fermions (i. e. at TV = 0, or A = 0) and 
when bi = 6 2 , if a is integrated out, A_ becomes con- 
strained to be zero. The model simply reduces then to 
3D HSE with two n-component complex fields. As dis- 
cussed in the previous section, this doubling of degrees 
of freedom should be considered to be an artifact of our 
relaxation of the y —> limit in the underlying lattice the- 
ory. We therefore impose that for N — our calculation 
reduces to the HSE with a single complex scalar field by 
formally taking the number of complex components for 
each field in our model to be n — 1/2, not n = 1. As 
we will be performing perturbative calculations, n will 
determine only the coefficients in our expansions, and 
therefore for n = 1/2 and N = those will be then ex- 
actly the same as in the HSE. This computational device 
will allow us to asses the effect of the coupling to gapless 
fermions. 

If we were to set 61 = 62 in the Eq. (6) we would 
have a global U{2n) symmetry associated with the 2n 
$ fields, when A_ = 0. However, the coupling to the 
gauge field A_ breaks such a U{2n) invariance. Hence, 
if initially bi = 62 , quantum corrections would cause bi 
and 62 to be no longer equal in the effective action. To 
consistently renormalize the theory we therefore need to 
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FIG. 1: Diagram showing non- vanishing terms that con- 
tribute to the effective theory obtained by integrating out 
the fermion fields. The smooth solid lines represent fermion 
propagators. The wavy lines represent gauge fields. 



keep track of both couplings in the intermediate stages 
of the calculation. Only at the very end will we impose 
the constraint that = ($2) in the effective action, 

once the quantum corrections have been computed. 



IV. RENORMALIZATION GROUP IN 3D 

We proceed to study the continuum theory in Eq. (6) 
by using one-loop renormalization group at fixed dimen- 
sion. Although such a calculation is essentially uncon- 
trolled, it has been previously used successfully in |l3j by 
appropriately choosing the renormalization conditions at 
the critical point t = 0. We will find that by increasing 
A one progressively destabilizes the second order phase 
transition in favour of the first order one. Above certain 
A = A c the phase transition becomes first order for all 
values of the couplings e, b\ and 62- To show that this 
may not simply be an artifact of our approximation, in 
the next section we reach a similar conclusion by using 
the minimal subtraction scheme in 4 — e dimensions. 

We begin by integrating out the fermionic fields and 
expand out in powers of the vortex gauge field a. Such 
an expansion is shown in Fig. I. For A <C 1 we may retain 
only the first non- vanishing term in such an expansion, 
which is proportional to the polarization tensor, IL^. 
The a field may then be integrated out exactly, and we 
are left with an effective theory 



s eff = 



d 3 a 



(V x A_(a;)) s 
2IVIA 2 



(8) 



where A 2 



is an effective "charge" associated 



with the A field, introduced by fermions. In a less ap- 
proximate treatment one would also find the terms of 
higher order in the gauge-field A_, as well as a renor- 
malization of the coupling in the A 2 , term. The effect of 
including these higher order corrections, will be to renor- 
malize the effective charge A, so that A will be a power 
series in both A and 1/N, with the value of A given above 
as the leading order term. For the sake of simplicity, we 
will however neglect these higher order corrrections. 

Next, we compute the polarization tensors 2nA 2 x~ 1/ 
and 2ne 2 x+ i ,, for the A_, A + fields, respectively, to one 
loop order, ft is easy to show that Xav = Xav to all orders 
in perturbation theory. We give an explicit expression 
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FIG. 2: Diagrammatic contributions to and Y ' . The 
smooth solid lines with numbers in them represent <&-field 
propagators, the number 1 corresponding to 3>i-field propa- 
gators, etc. The broken lines represent the propagators for 
A_, and the wavy lines for A+. 





FIG. 3: Terms that do not contribute to the RG flow for T^ 4 ] 

and The cross hatched circle represents any allowed 

combination <I>, A_, and A+ propagators. 



for in Appendix A. Another important observation 
is that there is no mixing between A_ and A+, again 
to all orders of perturbation theory. This is because A_ 
couples to $2 with opposite charge to the coupling of 
A to <&i, so contributions from $1 fields to the mixing 
term always cancel those of $2 fields. Finally, the 4>i 
propagator is exactly equal to the $2 propagator. In 
Appendix A we provide the expressions for the inverse 
3? propagator T^ 2 ). Contributions to x M „ and T^ 2 ) are 
shown diagrammatically in Fig. 2. 

We then consider the four point vertex functions for 
the 3? fields. For renormalization group purposes, we do 
not need to include all the terms that appear to one loop 
order. In the customary Landau (transverse) gauge, we 
adopt here, terms with the graphs of the form given in 
Fig. 3, do not affect the critical behaviour 0], since they 
vanish for finite t when p — > 0. We shall denote the four 
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point vertices for 3>i — $1 scattering and 3>2 



■* 2 byrW 



and r^, respectively. We shall use T^, to denote the 
four point vertex of scattering between two 4>i and two 
$2 fields. 

We calculate these vertices at the symmetric point 
Pi Pj =p 2 (4Sij - 1), 

where pi are the momenta of the external legs. Some 
care must be taken when working in any fixed dimension 
other than 4. The reason for this is that, when consid- 
ering T^l at the symmetric point, two of the scattering 
channels (the t- and u-channels (2^]) that contribute to 
the vertex have a total external momenta p\ — p2, instead 
of pi +P2 for the other (s-) channel. This implies, at the 
symmetric point, (pi +P2) 2 gets replaced by p 2 , as usual; 
but (pi —P2) 2 should be replaced by 2p 2 . So on calculat- 
ing the loop integrals, an extra factor of 2^ 2 ~ d / 2 ^ appears 
in the terms that come from the t- and u-channels. The 
same consideration must be taken also when calculating 
the four point functions I{ 2 an d ^21- Expressions for 



(4) 



r^, T 2 4 2 and at the symmetric point are given in 
Appendix A. For the purposes of illustration we show 
these terms in Fig. 4. 

To derive expressions for the renormalized couplings, 
we now impose the following renormalization conditions 
in 3D 

rW.*(p) = o, ^M = i, 

r M" R (Pl)-P2,P3,P4))sP = V bl,R, 
rf^V,P2,P3,P4))sP = ^, 



.(4) 



q=p/c 
p/c ' 



r( 3 )+^( P ,o,p)=p 1 /2 ei?5 

T(V-< R (p,0,p) = \ R . 



(9) 



Here, T^ + ' R and r^~< R are the renormalized A + — 
4>i — $1 and A — #1 — $1 vertices, respectively. As 
bare quantities, and r^ 3 ^ + take the values A and e. 

Because of the Ward identities (given in Appendix A), 
we do not require an expression for r^ 3 ^ + nor for r( 3 )-, 
to one loop order, to determine renormalized charges; all 
that is required are the expressions for -D+„ and D~ . 
Here, c is an arbitrary scaling factor which we shall tune 
|jl3f so that the value of the Ginzburg-Landau parame- 
ter bi /4e 2 at the tricritical point in standard HSE (i. e. 
when A = 0) matches the known value. Using the renor- 
malization conditions and the expressions in Appendix A 
we are able to write down the renormalized couplings for 
71= 1/2, d = 3 



pb\,R = b\- 



7b 2 , b\ 



32V2 16 32\/2 




r"= 




b.) 



FIG. 4: Terms that contribute to renormalized and Tj 4 ^. 



2 p 



e 2 b 2 



2A 2 6i 4e 2 A 2 



4\/2 P 
pbi,R = b 2 - 



16 



2p 4V2 P 



V2tt 2 
1 

+ 71 

4e 2 A 2 
V2tt 2 



4A 4 p 
V2tt 2 



16-\/2 



hb 2 



Pe 2 R 



4A 4 p 

V2^ 2 ' 

4 

e c 
~ 16p' 

~ ~16~ 

(10) 



From the last line one can see that the coupling A does not 
flow. In fact, the only reason why we found the A 4 term in 
Xr is that we are working right at the critical point t = 0; 
if this were not so the polarization would be proportional 
to p 2 and the coupling A would not renormalize at all. 
This is because the Maxwell-like term for A in the Eq. 
(8) is non-analytic in small momentum, and therefore can 
not get renormalized pfJ. We will therefore neglect the 
A 4 term in Xr, and consider A to be an exactly marginal 
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coupling hereafter. 

In the Eq. (10) we introduced e' as an infinitesimal 
quantity used to regularize the logarithmically singular 
terms that appear in the renormalized couplings, e' is de- 
fined to be an infinitesimal shift in the scaling dimension 
of the A_ gauge field propagator, so in 3D 



p 2 



We choose to regularize this way, as opposed to dimen- 
sional regularization, due to fact that, when we generalise 
our model to arbitrary dimension (in the next section), 
we find that these logarithmic poles are present in all 
dimensions. The reader can easily convince him/herself 
that by setting A = and b\ = b 2 = b one indeed finds the 
correct expressions for renormalized couplings for HSE 
to one loop order in 3D. On rescaling, these expressions 
agree with [JDj), except for one difference; the e 4 term that 
contributes to b\ and b 2 , contains an extra piece, missed 
in In 4D this extra contribution can be neglected, 
as it is not singular; but in 3D it gives al/p contribution, 
so must be included in the e 4 term. 

Before proceeding to the flow equations one needs to 
address the following technical issue: how does one treat 
the A 4 term in and &2,_r? If we were to blindly differ- 
entiate the expressions in Eq.(10), we would find the A 4 
terms in both the beta functions for the b\ and 62 fields. 
But, power counting in 3D tells us that such terms are 
regular (i. e. not singular in the infrared limit) and 
should be discarded (together with all other regular con- 
tributions) from the beta functions. On the other hand, 
as we shall see, when close to four dimensions these terms 
should be included as they are as relevant as all the other 
terms. We think that the correct way to handle such 
terms when working directly in 3D, is to absorb these fi- 
nite terms into bare values of the quartic couplings. That 
way they determine only the initial values of the flowing 
couplings (actually reduce b\ and 62) instead of entering 
the beta functions. At the end of this section we will 
briefly discuss what the effect of including such terms in 
the flow equations would be, to reassure the reader that 
our general conclusions would remain unchanged. 

By making such a finite shift of b\ and 62, and then 
differentiating the remaining expressions we arrive at the 
flow equations, 

dh __ 
^ dp 
,0 I 



bf 



7 

32^2 



1 
16 



32 V2. 
dh = 
dp 



16 



1 

71 



16^2 



bib 2 



-e 2 h 
2 


3 4 
4V2 


2A 2 6i 


4e 2 A 2 \ 


IT 2 


V2ir 2 ) ' 


-e 2 b 2 
2 


f % 4 + 
4^2 


2A 2 6 2 


4e 2 A 2 \ 




+ V2^ 2 J ' 


de 2 


4 

9 ce 


dp 


E + 16' 



A 


(&_,<*_ ) 


(b-,S+) 


(b+,8-) 


(64-,*+) 





(0.324,0) 


(3.03,6.44) 


(6.39,-5.53) 


(5.78,0) 


0.4 


(0.347,-0.060) 


(3.17,6.48) 


(6.51,-5.78) 


(5.83,0.289) 


1.0 


(0.455,-0.339) 


(3.90,6.67) 


(7.20,-6.92) 


(6.12,1.63) 


1.5 


(0.582,-0.660) 


(5.00,6.86) 


(8.24,-8.41) 


(6.55,3.40) 


2 


(0.711,-0.985) 


(6.64,6.81) 


(9.74,-10.3) 


(7.04,5.91) 


2.5 


(0.824,-1.27) 


N.A 


(11.7,-12.7) 


N.A 



TABLE I: Fixed points for various values of A. 



where we omitted the subscript R. 

First, we constrain the parameter c in the same manner 
as in |f|: we set A = and b\ — 62, calculate the posi- 
tion of the tricritical line in the flow diagram for the stan- 
dard HSE, and then use the lattice value of the tricritical 
Ginzburg-Landau parameter k c = 6i/(4e 2 ) = 0.42/V2 
[§7| to fix c. This yields a value of c=17.4. 

To analyse the flow equations it is convenient to define 
the new couplings b = bi, 5 = 62 — b\. In Fig. 5 we show 
the numerically calculated flow for e 2 = e* 2 = 16/17.4, 
the fixed point value for e. We see in Fig. 5, for A = 
and A = 2, that there are four fixed points. At A = we 
have the usual HSE fixed points at 5 = (b| . These are 
the tricritical point (6_ , (5_ ) , which lies on the k c line and 
has two repulsive eigenvalues in the e 2 = e* 2 plane; and 
the stable (critical) fixed point, (b + ,S + ), which has only 
attractive eigenvalues, and should be identified with the 
inverted XY critical point in our approximation. The two 
extra fixed points are for 5 0. Both have one attractive 
and one repulsive eigenvalue. The one for which 5 is 
positive we denote by (&_, <5+) and the one for which S is 
negative by (6 + , <5_ ) . By setting the beta functions equal 
to zero and solving the resulting equations one is able to 
find the fixed points for various values of A as shown in 
Table 1. Indeed, by linearizing our RG equations around 
these fixed points one is able to calculate the eigenvalues 
and eigenvectors and show that these coincide with the 
numerical solutions shown in Fig. 5. Looking at the 
Table 1 and the Fig. 5, we observe two important trends. 
For A 7^ 0, the fixed point (b + , 5+) is no longer at S = 0, 
but at a positive value of S. As A increases so does this 
value of 8. We also see that (6_,(S+) starts to move to 
larger values of b. At the value A c = 2.04, these two 
points cohere; and above this they move into the complex 
plane. For A = 1, this corresponds to N c 3.4 . We 
see in Fig. 5c that the effect of losing these two fixed 
points is to completely destabilize the second order phase 
transition. We are also able to show numerically that for 
the runaway trajectories, the coupling that appears in 
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FIG. 5: Flow diagrams in the (6, S) plane with initial condition 
e = e* for various values of A. The arrows show the direction 
of flow as p -> 0. a.) A = 0. b.) A = 2. c.) A = 2.5. 



the effective potential with (<I?i) = ($2), b + 5/2 — > 
—00 for p — > 0, which should lead to a fluctuation induced 
first order phase transition. 

In Fig. 6, we plot the separatrices in the space of 
(6, e 2 ) between the first and second order phase transi- 
tions for selected values of A. Notice that as A increases 
the boundary between the two phase transitions moves to 
the right, reducing the domain of attraction of the stable 
fixed point. 

If one would include the A 4 terms directly into the RG 
equations, the effect is to more rapidly destabilize the sec- 
ond order phase transition in favour of a first order tran- 
sition, by pushing down A c to A c = 1.14 (corresponding 
to N c — 1.05 for A = 1). Again the same flow structure 
is present, and (&+, 6+) is destroyed in the same manner 
as discussed above at A c . 

We next compute the anomalous dimension 77$ and the 
correlation length exponent v at the stable fixed point, 
for A < A c . These may be obtained by considering Z^ and 
Z^2 ; where the former can be calculated from the com- 

(2 1) 

posite Greens function T\ [ , for two separate $ fields 
and a <fr 2 insertion. The terms that contribute to one 
loop order are shown in Fig. 7. We find 



m = 



-,*2 



A 2 

T 2 



32 



4 7T^ 

e* 2 A 2 

4 7T 2 



(12) 



For N — 2 and A = 1 this implies v = 0.65 and 77$ = 
—0.48. Numerical values of the exponents at other values 
of A are provided in Table. 2. For increasing A, v slightly 
reduces in value and 77 becomes more negative. Notice 
that the value of v calculated for A = is very close to 
the known XY value of v = 0.667 & 




FIG. 6: The boundary in the e 2 — b plane between the fluc- 
tuation induced first order phase transition and the second 
order transition; calculated for A = 0.4, A = 0.8 and A = 1.2 
The straight line for A = is the k = k c line in the standard 
HSE H. 
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FIG. 7: Terms in F\ [ ' that contribute to Z^i to one loop 
order. 



A 


V 


77* 





0.663 


-0.230 


0.4 


0.663 


-0.246 


1.0 


0.659 


-0.331 


1.5 


0.652 


-0.458 


2 


0.637 


-0.635 



TABLE II: Critical exponents v and 77* calculated for various 
values of A. 



V. RENORMALIZATION GROUP NEAR 4D 

When we generalize our theory to arbitary dimension 
we need to make one modification to the Lagrangian in 
the Eq. 6, which is to replace the coupling A 2 by eA 2 . 
This is to insure that the theory has a well defined ex- 
pansion for small e = d — 4, and handles the effect of 
fermions in a controlled way close to 4D. On integration 
of the fermion fields we use the well known formula for 

n 

to 4D 



M „ for N massless fermion fields, to one loop order, close 



2 ,Np~ e 

IV = {p^Pv - d^p ) 12 ^ 2e 



0(1), 



to write down a generalised effective theory 
"(V x A_(.t)) 2 



S eff = 



d 3 x 



2eA 2 |V|( 4 - d ) 



(13) 



(14) 



where close to 4D, A 2 = A 2 N/12. From the expressions 
for the loop integrals in x^, T^ 2 \ r^ 4 j, and T^l 
given in Appendix A, for arbitrary dimension, it is easy 
to extract terms singular in e = 4 — d. It should be 
pointed out that as well as including e-poles we also have 
to include e'-poles (e' is defined in the previous section) 
coming from terms that contain logarithmic singularities 
in all dimensions, due to the form of the effective A_ 
propagator. Since we are using the minimal subtraction 
scheme, we define Z$, Z^, Za+, Za-, and Z^ so 
as to subtract out the only e and e' poles in the loop 
integrals. We get the following forms for the renormalizcd 
couplings, 



3e 2 p~' 



&f(n + 4) bin 



V 



+ 



3e\ 2 p- 
4ttV 

6eAV 



bi 



87r 2 e 

6e 4 p _e 
47r 2 e 



4tt 2 



+ 



6e 2 eA 2 p 
An 2 £ 



9 



b 2 ,RP e = b 2- [J>2 + 6i6a(n + l)] 



P 



3e 2 p~ ( 
47r 2 e 



+ 



3e\ 2 p-' 
47r 2 e' 

6eAV 
4tt 2 



87r 2 e 
6e p~ e 
47r 2 e 



6e 2 eA 2 p 
47r 2 e 



4 - 

ne p 
~L2^ 



(15) 



From the above expressions we derive the following beta 
functions 



dbi , 
p-j— = -eoi + 
dp 



I 



&5™ 



1 

8^ 



3(e 2 + eA 2 )6i , 6(e 4 + e 2 eA 2 + e 2 A 4 ) 
4tt 2 
db 2 

p— = -eb 2 
dp 

3(e 2 + eA 2 )6 2 6(e 4 - e 2 A 2 + e 2 A 4 ) 



4tt 2 

b\ + hb 2 f y + 1 



4tt 2 



4tt 2 



P- 



dp 



I 4 

n e 
24^ : 



(16) 



where n' — 2n is the total number of complex fields in 
3>i and 4>2- For the two fixed points (eb*, eS*, ee* 2 ), that 
become the charged fixed points in the standard HSE for 
S = 0, A = 0, we find the following equations for small A 



16^ 2 



(2n' + 8) 



1 



3(e 



*2 



A 2 



S*n' 



H 

3(2n' 



4tt 2 16tt 2 
3(e* 2 + A 2 ) S*n' \ 



16tt 2 
3e* 2 A 2 



4tt 2 

(e* 4 + e* 2 A 2 ) 



167T 2 J 

1/2 



(17) 



where e* 2 = . The fixed points are real for n! > n' c , 
when it is possible to have a second order phase transition 
in the theory. n' c is therefore determined by 



3(e* 2 + A 2 ) S*n' c \ 2 3(2< + 8) 



4tt 2 



16tt 2 



16tt 2 



(e**+e* a A a ) = 
(18) 

as well as with the Eq. (17). We solve these equations 
to second order in small A: 



n' c = 182.9 + 126.9A 2 + 0(A 4 ). 



(19) 



For A = we recover the well known result of JL4| . When 
small A is present n' c increases. This we take as an in- 
dication that again the second order phase transition is 
destabilized by fermions, in qualitative agreement with 
the conclusion of the previous section. 



VI. NEUTRAL THEORY IN 3D 

It is interesting to see what the coupling to fermions 
does in the "neutral" case e = 0, which would reduce 
to the XY (not to the inverted XY) model when A = 0. 
Such a theory would correspond to the charged XY model 
for the fluctuations of the phase of the superconducting 
order parameter, which is dual to the Higgs scalar electro- 
dynamics with a massive gauge field A + |E3[. Assuming 
that such a gauge-field mass is large one can neglect the 
A+ field altogether in the Eq. (6), and define what we 
will call the "neutral" version of the theory. 

Setting A + = in the Eq. (6) allows one to greatly 
simplify the theory by making the transformation & 2 — > 
*&2; upon which the theory becomes symmetric under 

4>i < ► $2- One may then set 3>i = $2 = * and 

bi = b 2 = b from the outset, and take $ to have one 
complex component. This leads to £ = C' h + £/ + £ c , 
with 



4 = |(V-zA_)*| 2 + t |*| 



1*1 



(20) 



A closely related model has been discussed before in [^9| . 
The beta function for the single coupling constant b is 
readily obtained by taking the e = limit of the Eq. 
(11): 



P 



db 
dp 



-b- Aa 2 6+(2V2 + 1)^. (21) 

TT Z 16 



In contrast to Eqs. (11), now there is always a stable 
fixed point, at which 



A 2 

m = -— , 

7T Z 



^2 ■ 



(22) 



(23) 



For A = 1 and N = 2 (A = tt 2 /4) this yields 77$ = -0.25 
and v = 0.68. 

There is still, however, a possibility of a discontinu- 
ous transition in the theory, since upon integration over 
A_ the bare value of coupling b becomes reduced by a 
contribution ~ A 4 , similarly as in the Eq. (10). If b is 
small this will turn it negative and the transition should 
become discontinuous. For large enough b, on the other 
hand, the flow is towards the modified XY critical point 
with the above exponents. 

A related theory of coupled Dirac fermions and scalar 
field that also exhibits iV-dependent exponents has re- 
cently been considered in [M. 



VII. CONCLUSION AND DISCUSSION 

To summarize, we have argued that the coupling of vor- 
tex loops to gapless fermionic excitations tends to desta- 
bilize the continuous inverted XY transition in favour of 
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a fluctuation induced first-order one. With increase of 
the coupling to fermions, the domain of attraction of the 
inverted XY critical point in 3D scalar electrodynamics 
shrinks, until it finally disappears at some critical value. 
In our approximation the critical number of fcrmion com- 
ponents where this happens is N c as 3.4, which suggests 
that for the physical case N = 2 the quantum phase 
transition between dSC and the SDW may still be con- 
tinuous. We showed that the correlation length exponent 
decreases with the number of fermions and estimated it 
to be v rj 0.65 for the physical case N = 2. 

Our number for N c , which happens to be larger than 
two in our approximation, may well be an overestimate. 
There are two further effects which can reduce it. First, 
as we already discussed, retaining A 4 terms in the beta 
function would yield a lower N c . More importantly, 
tuning our free parameter "c" in the beta functions for 
HSE to a larger value of the tricritical Ginzburg-Landau 
parameter would also decrease N c . For example, for 

k c = 0.8/V2 @, H, one nnds that N <= is reduced all 
the way to zero, and the transition becomes discontinu- 
ous for all values of the coupling constants. 

While there is a great degree of uncertainty in the pre- 
cise value of the critical number of components N c inher- 
ent to our calculation, the qualitative result that coupling 
to gapless fermions favours a discontinuous transition, 
we believe, is likely to be correct. This is because this 
coupling introduces another massless gauge field into the 
HSE, the effect of which is, crudely, to make the dual 
theory effectively more type-I. The reader should recall 
that, in the phase diagram of standard Higgs scalar elec- 
trodynamics jyj), increasing the coupling to the gauge- 
field and making the system more type-I would always 
eventually lead to the first-order transition. The effect 
of the additional gauge field is to increase the type-I re- 
gion in the space of coupling constants at the expense 
of the type-II, where the transition is continuous. The 
lowest order calculation suggests that above certain crit- 
ical value of this coupling, that may be measured by the 
number of fermionic components N, the type-II region 
disappears completely. 

To put our calculation in perspective, it may be useful 
to ask what the result would be if the superconductor 
were not gapless, but had a full gap. One can model 
that situation by adding a mass term ~ m^fj^fj to £/. 
Integration over fermions would then give an ordinary 
Maxwell terms for a, which when integrated out would 
finally yield a massive term for A_, mA. 2 _ t , where by 
subscript t we denote the transverse components. Such a 
massive gauge-field is irrelevant, and the transition would 
remain in the inverted XY universality class. (Similar re- 
sult is found in the gapless case but at T 7^ 0, where finite 
temperature provides an effective mass for the gauge field 
p| , and makes it irrelevant at the standard KT critical 
point.) We find it pleasing that in this case one recovers 
the result that intuitively appears to be correct. 

In the previous work JtJ the procedure we followed was 
complementary to the one in this paper: there we inte- 



grated approximately the vortex degrees of freedom, to 
find that when vortices are condensed electrons organise 
themselves into a SDW. This, however, left the question 
of the nature of the dSC-SDW transition open. Here 
by integrating fermions first we study the nature of the 
phase transition from the vortex point of view. We em- 
phasise, however, that one is still considering one and the 
same phase transition in the theory (1). We expect the 
T = free energy defined by Eq. (1) to have only one 
singularity as the parameter K is varied, at which vortex 
loops blow up (and superconductivity is lost), and simul- 
taneously chiral symmetry for fermions becomes broken 
(and the SDW is established). This is precisely what 
happens in the mean-field treatment of the vortex fields, 
and we believe should be generally true. Our expectation 
is apparently met in the lattice calculations on a related 
model of fermions coupled to a bosonic fields via gauge 
fields |33| , where indeed only a single phase transition is 
found ]p4| . 

Theory similar to our Eq. (6) has also been recently 
studied by Ye |55|, but with one important difference: 
there is only a single scalar field while $2 = 0. In 
such a theory A + and A_ are coupled to each other, 
and in particular, when =/= the integration over 

A + ultimately makes a massive in the insulating phase. 
As discussed elsewhere 0, this is an artifact of the par- 
ticular singular gauge transformation made in arriving at 
the theory. Here we only note that by demanding that 
($1) = ($2), as dictated by the underlying lattice model, 
completely decouples the gauge fields A + and A_, which 
then leads to the massless gauge field in the insulating 
phase as in Eq. (7). 

It may be interesting to note a parallel between our 
model and another popular field theory, namely the 
Thirring model |3(| . Assume that one in the super- 
conducting state, to > 0, and integrate all the fields in 
(6) other than fermions, in the Gaussian approximation. 
Since in the superconducting state the gauge field a is 
massive, one finds a quartic term 

in the remaining action for fermions. The theory with 
just this quartic term would be the 3D Thirring model 
with N = 2. This model has a transition into a phase 
with broken chiral symmetry as to is decreased. While 
our field theory is not exactly equivalent to Thirring 
model |37j] , it is interesting that we find the same trends 
as observed in its numerical studies p^| : a) the correla- 
tion length exponent v decreases with N, b) above cer- 
tain N c the phase transition becomes discontinuous. It 
would be interesting to see if these non-trivial features of 
the Thirring model can be obtained in some RG scheme 
similar to ours. 

As applied to underdoped cuprates, our calculation 
would suggest that the T 7^ line of the Kosterlitz- 
Thouless transitions ends in a modified XY critical point 
with new critical exponents, or even in the first-order 
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transition, at T = 0. The latter would imply, for ex- 
ample, that T = supcrfluid density p s f{T = 0) may 
discontinuously drop to zero as doping is reduced, while 
the former would mean that p s f(T = 0) <~ (x — x c ) z ", 
where x c the critical doping, and z = 1 the dynamical 
critical exponent. Detection of these effects is at present 
probably prohibited due to disorder that is most difficult 
to control precisely in the underdoped region. In the 
mean time, numerical calculations on the lattice model 
(1) should be able to shed some additional light on the 
problem. 
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APPENDIX A: EXPRESSIONS FOR VARIOUS 
1PI VERTICES TO ONE LOOP ORDER 

What follows is a general derivation of various 1PI ver- 
tices to one loop order in arbitrary number of dimensions. 
By applying Feynman rules to the effective theory (de- 
scribed by Eq. (14) in the text) we find the following 
forms for the polarization tensor 



d d k 



(27T) 



~ d G{k) 



d d k 
(2^ 



(fc„ + (k + P ) tl )G{k)G{k - P )(K + (k+ p) v ) 



(Al) 



Here, G(p) denotes the bare greens function for both $i 
and * 2 fields. For we find the following 



r?\p) = G( P y 



+2e 



d d k 



&i(n + l) b 2 n 
2 + ~Y 



d d k 

W) 



-Mk) 



(27r) d 



(k) 



(2ir) d 

/ i0yi( k +P)» D tA k -p)(k+p)vG{k) 



d d k 
J2^ a 



(k+ P )^D+Jk-p)(k+p) v G(k), 



(A2) 



where D+ (k) and D~ (k) are the bare A_ and A + 
gauge field propagators, respectively. For T^\, Tfl and 
we find 



r (4) 
1 1,1 



.(4) 
2,2 



b 2 (n + 3) + nb\ b\ 

23-d/2 + 2 



d d k 
(2tt)° 



G{k)G{k + p) 



_ 2 d/2 e 4 



-2 d ' 2 e 2 e\ 2 



_ 2 d/2 e 2~ x 4 



r (4) _ h 

L l,2~ 2 



d d k 

J2^Y 
d d k 



D+„(k+p)D+ v (k) 



—D+Ak+p)D-Jk) 



d d k 



—D- v (k + p)D- v (k) 
bib 2 (n + l) 



bl b 2 

— A h 

4 x 2 3 ~ d / 2 



2 2-d/2 



I 



d d k 

j^y d G(k)G(k +P ) 



-2 d ^J ^0- d Di v (k+p)Dl u {k) 
+2 d / 2 e 2 e 2 l 2 j 0^Dl v (k+p)D- v (k) 
- 2d ' 2 ^j jd D- A k+p)D-Ak). 



(A3) 



We start by evaluating the following by dimensional 
reglarization 



d d k 
(27T) d 



(2tt) c 



(fc M + (k+p)Mk)G(k-p)(k v + {k+p) v ) 



< f ddk ( — 

7 J2^y \¥ 



{2k„ + p„){2k v + p v ) 



^k 2 k 2 {k — p) 2 

2r(i-d/2)r(rf/2) 2 , 4 , 2 



T{d){An) d / 2 



P (p 5^-p^pv), 
d d k 1 1 



d d k 
(2^p 



y k 2 (k + p) 2 

T{2-d/2)T{d/2-l) 2 d _ A 
T(d- 2) (47^/2 P ' 

{k+p)^D+ (k-p)(k+p) v G(k) 



d d k {k 2 p 2 -(k-pf 



= (d-l) 



(27r) d k A {k+p) 2 
p" 



T(2-d/2)T(d/2-l) 2 . 4 



T(d- 2) (4ir) d / 2 
d d k 



I ^dDUk+P)DtM^ 



d d k 
(2^p 



5uu 



k^k v 



(k+p) fl (k+p) L 

(k+p) 2 



d T(2-d/2)T(d/2-l) 2 d _ 4 

4 l ' T(d-2)(4n) d / 2 P • 



(A4) 

Terms involving a single D~ v propagator in both x^v 
and the vertices are logarithmically divergent for all 
numbers of dimensions, so as discussed in in the text we 
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regularize these terms in by introducing a shift e' in the 
scaling dimension of the A__ propagator. Using this, we 
write these integrals in the following way 



+ p)„D- v {k - P )(k+ p) u G{k) 

d d k {k 2 P 2 - (k. P ) 2 )k 4 - d - e ' 
(2ir) d k 4 (k+p) 2 

2(d-l) p- e ' 



T{d/2){A-K) d / 2 e' ' 



f & k u k u \ f & {k + p)^(k+p) v 



(2n) d \ ^ k 2 7 V (k+p) 2 

k 4 ~ d ~ e ' 2(d-l) p 



k 2 (k+p) 2 {4n) d / 2 T(d/2) e> 
We also evaluate the A 4 term in both 61 and 62 



(A5) 



(27T 



D-J^D-Jh+p) 



{5 k„k v \ f & {k + p)^{k+p) u 



(27r) d V k 2 J \ ttv (* + P) 2 

k 4 ' d \k+p\ i - d 
k 2 {k+p) 2 

(d-l) ( T(d/2-2) r(d/2-l) \ ,_ d 
(47r) rf / 2 \r(d/2- l) 2 2T(d/2) 2 J P ' 



(A6) 



By using the RG conditions given in the text (Eq. (9)) 
and the Ward identities 



(k - p)^ + ' R (k, k~p,p)= e(T^ (k) - T ( V (p)), 

(k -p)„rW-> R (k, k-p, P ) = \(r( 2 \k) - r^( P )), 



(A7) 



one arrives at the Eq. (10) for the renormalized couplings 
in 3D and for n = 1/2. 

Close to 4D and at the critical point, we have the fol- 
lowing singularities 



Lw(p)L 





-A 




sing \ 



1 p- 

24tt 2 e 
3e 2 p~ 



■(P 2 5y.v -PfiPv), 

3eA 2 p- e ' \ 



8ir 2 



t?Up) 



r 6 2 (n + 4) bin 



sing 



8vr 2 6 



6e 4 p e 
47r 2 e 

= b_2 

sing 2 

6eV 



6eA 4 p 



4tt 2 



6eA e l p 
47r 2 e' 



6| 



2 

6eA 4 £T 



6x62(72+ 1) 



87T 2 6 



6eA e p 



47T 2 6 



47T 2 47T 2 e' 

From these it is straightforward to define Z$, Z^i, Za+, 
Za~, Zbi and Z(,2, that remove these singular terms. 
Using the relations between renormalized and bare cou- 
plings e\ = p- e Z A+ e 2 , \ 2 R = Z A -\ 2 bin = p-'Z 2 Z^bi 
and &2_r = p~ e Z^ > Z^ 2 1 b2 we arrive at the Eqs. (15) given 
in the text. 
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